Zoological Research 35 (5): 362-372 DOI: 10.13918/j.155n.2095-8137.2014.5.362 





Influence of a large dam on the longitudinal patterns of fish assem- 
blages in Qingyi Stream 


Xiao-Yun SUI'?, Zhi LU’, Yun-Zhi YAN, Yi-Feng CHEN", Yin-Tao JIA? 


1. School of Life Sciences, Peking University, Beijing 100871, China 
2. Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan 430072, China 
3. College of Life Sciences, Anhui Normal University, Wuhu 241000, China 


Abstract: Using seasonally collected data (2009—2010) from 15 sampling sites that represent first- to fifth-order streams within 
the Qingyi watershed, we examined the spatio-temporal patterns of fish assemblages along two longitudinal gradients to explore 
the effects of a large dam on fish assemblages at the watershed scale. No significant variation was observed in either species 
richness or assemblage structure across seasons. Species richness significantly varied according to stream order and gradient. 
Dam construction appeared to decrease species richness upstream substantially, while a significant decrease between gradients 
only occurred within fourth-order streams. Along the gradient without the large dam, fish assemblage structures presented 
distinct separation between two neighboring stream orders, with the exception of fourth-order versus fifth-order streams. 
However, the gradient disrupted by a large dam displayed the opposite pattern in the spatial variation of fish assemblages 
related with stream orders. Significant between-gradient differences in fish assemblage structures were only observed within 
fourth-order streams. Species distributions were determined by local habitat environmental factors, including elevation, substrate, 
water depth, current discharge, wetted width, and conductivity. Our results suggested that dam construction might alter the 
longitudinal pattern in fish species richness and assemblage structure in Qingyi Stream, despite the localized nature of the 


ecological effect of dams. 
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Stream fish assemblages are structured by abiotic 
factors, biotic interactions, and historical processes 
(Hoeinghaus et al, 2007). The physicochemical stream 
environment exhibits spatial heterogeneity and temporal 
variability, and may cause spatiotemporal variations in 
fish assemblages (Meador & Matthews, 1992). Typically, 
among the multiple spatial factors affecting fish 
assemblages, longitudinal-gradient variations in environ- 
mental conditions substantially influence the distribution 
and abundance of stream fishes (Araújo et al, 2009; 
Inoue & Nunokawa, 2002; Taniguchi et al, 1998; Vann- 
ote et al, 1980). From headwaters to downstream, fish 
species richness generally increases due to increases in 
stream size (Matthews, 1986), habitat diversity (Gorman 
& Karr, 1978), and shelter availability (da Silva Abes & 
Agostinho, 2001), and the different rates of fish immig- 
ration and extinction (Power et al, 1988). Concurrently, 
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fish species composition may vary longitudinally by 
species addition and/or species replacement (Boys & 
Thoms, 2006; Gorman & Karr, 1978; Roberts & Hitt, 
2010). The River Continuum Concept (RCC) describes 
the changes in community structure and species richness 
of organisms from headwaters to mouth waters and 
relates these changes to flow regime, water temperature, 
food availability, and substrate conditions (Vannote et al, 
1980). However, due to various human activities (e.g., 
damming, pollution, erosion), few riverine ecosystems 
remain free-flowing over their entire course. Overlaying 
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this pattern displayed by the RCC, the Serial Discon- 
tinuity Concept (SDC) claims that regulation by dams 
disrupts the underlying continuum and produces a series 
of lentic and lotic reaches, thereafter causing longitudinal 
shifts in abiotic and biotic parameters and processes 
(Ward & Stanford, 1983). 

Dams impact stream fishes and invertebrates in div- 
erse ways, such as blocking migratory pathways (March 
et al, 2003), fragmentizing habitats (Travnichek et al, 
1993), altering natural flow regimes (Bonner & Wilde, 
2000) and food webs (Power et al, 1996), decreasing 
water temperature downstream (Clarkson & Childs, 2000) 
and current velocity upstream (Bennett et al, 2002), 
disrupting riparian plant communities (Nilsson et al, 
1997), and shifting water chemistry (Humborg et al, 
1997). The extent to which stream fishes and invert- 
ebrates are affected by dams may be associated with the 
characteristics of dams (e.g., location, purpose, and 
management) and fauna (Cumming, 2004; March et al, 
2003). For example, large dams without spillways are 
impermeable barriers for migratory organisms and may 
extirpate all native migratory fish from upstream habitat 
(Holmquist et al, 1998). In contrast, large dams with 
spillways provide a possible passage for some native 
fishes. To date, investigations have primarily identified 
the effects of dams on lotic reaches directly below dams, 
mainstream reservoirs directly above dams, and lotic 
reaches upstream of impoundments (Clarkson & Child, 
2000; Cumming, 2004; Santucci et al, 2005; Travnichek 
et al, 1993). However, little attention has been paid to the 
effects of dams at the watershed scale. 

Stream ecosystems in China are experiencing mas- 
sive ecological perturbation due to diverse human 
activities, such as extensive agricultural and industrial 
production, urbanization, development, and hydraulic 
engineering construction, all of which threaten freshw- 
ater fish (Chen, 2005). Based on data collected from 
first- to fifth-order streams, representing two longitu- 
dinal gradients (one disrupted by a large dam and the 
other not) within a tributary (Qingyi Stream) of the 
lower reaches of the Yangtze River, the spatial and tem- 
poral patterns of fish assemblages were examined. Tota- 
lly, our three key goals were to: (1) to determine how 
fish assemblage patterns vary longitudinally from first- to 
fifth-order streams; (2) to identify the correlations 
between local habitat and fish assemblages; and (3) to 
assess how a large dam affects the longitudinal patterns 
of fish assemblages at the watershed scale. 
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MATERIALS AND METHODS 


Study area 

Qingyi Stream originates from the northern portion 
of Huangshan Mountain, and flows northeast toward its 
confluence with the lower reaches of the Yangtze River. 
This watershed is approximately 309 km long and covers 
an area of 7195 km’. Due to the subtropical monsoon 
climate, temperature and precipitation in this area are 
quite asymmetric across seasons. Annual air temperature 
ranges from —2.1 °C (January) to 27.5 °C (July), with a 
mean of 17.8 °C. Annual precipitation is high (approxi- 
mating 2000 mm/year), with most rainfall (79%) occur- 
ring in spring and summer (April to September) and only 
minor rainfall (less than 5%) occurring in the cold, dry 
winter (December to February) (Yan et al, 2011). 

Chencun Hydropower Station was constructed at the 
mid-reaches of the Qingyi mainstream in the 1950s. It is 
the largest station along the lower reaches of the Yangtze 
River, with 119 m deep impoundment and covering an 
area of 88.6 km’. Derived from Anhui Province 
topographic maps (1:300 000 scale), the Qingyi Stream 
is categorized as fifth-order (Strahler, 1957) at its 
mainstream. A total of 75 first-order, 34 second-order, 12 
third-order, 4 fourth-order, and 1 fifth-order stream are 
included in this watershed. Among the fourth-order 
streams, the Qingxi, Shuxi, and Machuan streams flow 
into an artificial reservoir (1.e., Taiping Lake) formed by 
the impoundment of Chencun Hydropower Station, while 
Huishui Stream directly flows into the mainstream with 
its confluence downstream of the station (Figure 1). 


Fish sampling 

A total of 15 segments representing first- to fifth- 
order streams were surveyed seasonally during April, 
August, and October 2009, and January 2010. The 
sampling segments were derived from two longitudinal 
gradients. One gradient (A) was set from the headwaters 
of Qingxi and Shuxi streams to the lower reaches of the 
mainstream, representing the gradient disrupted by 
Chencun Hydropower Station, and the other gradient (B) 
was set from the headwaters of Huishui Stream to the 
mainstream, representing the gradient not disrupted by a 
station. Each surveyed segment was set on the 
topographic map based on their stream orders. Four 
segments ranging from first- to fourth-order streams 
were selected within Qingxi, Shuxi, and Huishui streams, 
respectively. Three other segments representing fifth-order 
streams were set at the mainstream (Figure 1). Sampling 
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Yangtze 






Anhui Province 


Figure 1 Map of the Qingyi Stream watershed in Anhui Province, China 


Solid circles mark the 15 sampling sites. Sites 1-3, 4—6, 7-9, 10-12, and 13-15 represent first-, second-, third, fourth-, and fifth-order streams, respectively. 


Sites 1, 2, and 3 are located at Qingxi (QS), Shuxi (SS), and Huishui Streams (HS), respectively. »-< shows Chencun Hydropower Station. 


sites were selected in the field based on habitat represe- 
ntativeness and accessibility. Each shallow-water (<1.0 m 
depth) site of first-, second-, and third-order streams was 
sampled using backpack electro-fishing by wading in two 
passes. Each site encompassed at least two mesohabitat 
units (pool and riffle). Each deep-water (71.0 m depth) site 
of fourth- and fifth-order streams was sampled twice using 
boat electro-fishing along both riversides. Sampling sites 
were far (721.0 km distance) from conspicuous human 
disturbance, such as dam, farmland, and urban land, and 
each sampling was conducted with comparable effort (i.e., 
100 m long and approximately 40 min). Fish were 
identified to species level (except for Ctenogobius due to 
the deficiency in identifying tools), counted, and returned 
to sampling sites if alive. Voucher specimens were placed 
in 8% formaldehyde solution for further identification. 


Environmental survey 
Each sampling site was characterized using 10 


variables to describe local habitat: elevation (m), wetted 
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width (m), water depth (m), water temperature (°C), pH, 
conductivity (mS/cm), dissolved oxygen (mg/L), current 
velocity (m/s), discharge (m/s), and substrate. Elevation 
was determined by a portable GPS receiver. Wetted 
width was measured along five transects, regularly 
spaced across the stream channel. Water depth, temperat- 
ure, dissolved oxygen, pH, and conductivity were 
surveyed at four equal-interval points along each transect. 
Water depths «1.0 m were measured using a graduated 
wading rod, and depths 71.0 m were measured using a 
supersonic echo sounder. Current velocity was taken at 
60% of water depth at each point of each transect. 
Discharge of each channel was determined at the transect 
that yielded the most accurate measurement (smooth 
bottom and laminar flow). Along each transect, the 
proportion of substrate categories (particle size 1= 0-1 mm; 
2=1-5 mm; 325-25 mm; 4-25—-50 mm; 5-50-100 mm; 
6-100—500 mm; 7-500-1 000 mm; 8=>1 000 mm) 
was visually estimated, and an index of substrate coars- 
eness ranging from 1 to 8 was derived for each site 
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following Bain et al (1985). 


Data analysis 

Each fish species collected in this study was 
examined for frequency of occurrence (F) and relative 
abundance (RA), estimated from F;-100(S/S)?o and 
P=100(n/n)%, respectively, where S; and S are the 
abundance of the samplings of which species i were 
collected and of total samplings, respectively, and n; and 
n are the individual numbers of species i and total fish 
species, respectively. The above F and P values were 
determined independently for three groups of sampling 
segments, i.e., first- to fourth-order streams along 
gradient A and B and fifth-order streams, respectively. 
Species richness and fish abundance were analyzed using 
a separate three-way ANOVA model with stream order 
(first- to fifth-order), season (January, April, August, 
October), and dam (gradient A and B) as factors. One- 
way ANOVA was used to test differences in species 
richness across stream orders along gradients A and B, 
respectively. The Newman-Keuls test was used for post- 
hoc comparisons after ANOVA. Independent sample t- 
tests were used to compare species richness within the 
same ordered streams between the two gradients. All data 
were log-transformed to meet the assumptions of 
normality and homogeneity of variances. 

Discrete spatial patterns in fish assemblages were 
identified using PRIMER 5 (Primer-E Ltd., 2001). 
Following a Bray-Curtis similarity matrix calculation, an 
analysis of similarity (ANOSIM) was used to test 
variations in fish assemblages across stream orders, 
seasons, and gradients. Firstly, variations associated with 
stream orders and seasons were analyzed separately for 
gradients A (two-way crossed ANOSIM) and B (one-way 
ANOSIM due to no replication), respectively. Secondly, 
variations associated with gradients and seasons were 
analyzed separately for each stream order using two-way 
crossed ANOSIM. The relationships among assemblages 
from each site and season were graphically presented 
using non-metric multi-dimensional scaling (NMS) 
analysis. The contribution of each species to differences 
among assemblage groups was identified using similarity 
percentages analysis (SIMPER) (Clarke & Warwick, 
2001). The variability in fish assemblages in relation to 
local habitat environment was evaluated by canonical 
correspondence analysis (CCA) using CANOCA 4.5. All 
variables entered the CCA after a forward-selection 
procedure, showing their importance in explaining total 
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variability in species composition. The significance 
(P<0.05) of the CCA gradients was assessed by Monte 
Carlo permutation tests, and their importance measured 
by the eigenvalues of the first two axes (ter Braak & 
Verdonschot, 1995). All variables of fish assemblages 
(i.e. species richness and abundance) and habitat 
log(X+1) transformed 
assumptions of multivariate normality and to moderate the 


environment were to meet 
influence of extreme data. Species that occurred in less 
than two sites were excluded from the above analysis to 


avoid negligible weighting (Gauch, 1982). 


RESULTS 


Overview of species diversity 

A total of 13 647 fish were captured throughout this 
study, representing 57 species, 15 families, and 5 orders. 
Species of family Cyprinidae comprised on average 
59.6% of total species richness Species richness per 








sampling site amounted to 11.9+6.6 (mean+SD) species, 





and abundance per site was 227.6+205.5 specimens. A 
total of 16, 24, 24, 47, and 44 species were collected in 
first- to fifth-order streams, respectively. Cyprinus carpio, 


Mylopharyngodon piceus, | Hemiculter  leucisculus, 
Sarcocheilichthys | sinensis, Abbottina obtusirostris, 
Macropodus chinensis, Mystus macropterus and 


Hyporhamphus intermedius were only collected within 
fifth-order streams. Within the first- to fourth-order 
streams, 18 species were only collected from streams 
along gradient B, including Distoechodon tumirostris, 
Elopichthys bambusa, and Spinibarbus hollandi. Sauro- 
gobio dabryi was a unique species occurring in gradient 
A but not gradient B. Zacco platypus, Acrossocheilus 
fasciatus, Pseudogobio vaillanti, Vanmanenia steno- 
soma and Ctenogobius sp. were common (>40% of F) 
and relatively abundant (>1% of P) within both gradi- 
ents. In addition, Acheilognathus barbatulus, A. chan- 
kaensis, Rhodeus ocellatus, Squalidus argentatus and 
Misgurnus anguillicaudatus were more frequent and 
abundant within gradient A, while Opsariichthys bidens, 
Hemiculter bleekeri, Sarcocheilichthys parvus, Parab- 
otia fasciata, Siniperca chuatsi, Odontobutis obs- 
curus, Silurus asotus and Pseudobagrus albomargin- 
atus were more frequent and abundant in gradient B 
(Table 1). 


Local species richness 
Fish species richness significantly varied according 
to stream order and gradient, but not season (Table 2). 
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Table 1 Frequency of occurrence (%; left of dash) and relative abundance (%; right of dash) of fish collected within Qingyi 








Stream 
First to fourth-order Fifth-order 
Species Code 
Gradient A Gradient B 
Cypriniformes 
Abbottina rivularis ABR 50.0/1.7 31.3/2.3 83.3/2.9 
Abbottina tafangensis ABT* 3.1/«0.1 6.3/0.1 
Abbottina obtusirostris ABO* 16.7/0.2 
Acheilognathus barbatulus ACB 53.1/3.9 12.5/0.3 58.3/0.8 
Acheilognathus chankaensis ACC 18.8/0.9 6.3/0.3 41.7/0.5 
Acrossocheilus fasciatus ACF 68.8/7.7 87.5/10.1 25.0/0.3 
Carassius auratus CAA 28.1/1.8 31.3/2.1 91.7/5.8 
Cobitis rarus COR 53.1/5.6 37.5/2.2 
Cobitis sinensis COS 28.1/0.6 18.8/1.3 8.3/«0.1 
Cyprinus carpio CYC* 16.7/0.1 
Culter erythropterus CUE* 6.3/«0.1 
Distoechodon tumirostris DIT* 6.3/0.1 
Elopichthys bambusa ELB* 6.3/«0.1 
Erythroculter ilishaeformis ERI 6.3/0.2 41.7/2.3 
Gnathopogon taeniellus GNT 18.8/0.6 25.0/0.4 
Gobiobotia tungi GOT* 12.5/0.1 
Hemiculter leucisculus HEL 75.0/7.9 
Hemiculter bleekeri HEB 6.3/0.2 31.3/4.3 58.3/1.2 
Hemibarbus labeo HEL 18.8/0.4 
Hemibarbus maculatus HEM 31.3/0.7 16.7/0.1 
Misgurnus anguillicaudatus MIA 37.5/1.1 18.8/0.4 8.3/«0.1 
Mylopharyngodon piceus MYP 41.7/0.8 
Opsarüchthys bidens OPB 21.9/1.0 56.3/1.6 66.7/2.1 
Phoxinus oxycephalus PHO 12.5/0.2 18.8/3.8 
Pseudobrama simoni PSS 18.8/1.6 58.3/3.8 
Rhodeus ocellatus RHO 65.6/13.9 31.3/2.3 66.7/3.8 
Spinibarbus hollandi SPH* 12.5/0.1 
Zacco platypus ZAP 100/34.7 93.8/26.0 58.3/5.5 
Parabramis pekinensis PAP 6.3/0.1 16.7/0.1 
Parabotia fasciata PAF 18.8/0.6 31.3/2.0 33.3/5.0 
Pseudorasbora parva PSP 15.6/0.6 25.0/2.7 50.0/2.2 
Pseudogobio vaillanti PSV 43.8/2.5 43.8/4.8 33.3/0.5 
Sarcocheilichthys sinensis SAS* 8.3/«0.1 
Sarcocheilichthys parvus SAP 3.1/0.1 25.0/2.9 66.7/2.8 
Sarcocheilichthys nigripinnis SAN 25.0/4.0 75.0/2.3 
Squalidus argentatus SQA 344/1.2 12.5/0.7 33.3/1.0 
Saurogobio dabryi SAD 6.3/0.1 8.3/«0.1 
Varicorhinus barbatulus VAB* 12.5/0.2 
Vanmanenia stenosoma VAS 65.6/5.2 50.0/3.9 25.0/1.1 
Perciformes 
Channa argus CHA 18.8/1.4 91.7/6.1 
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(Continued) 
First to fourth-order Fifth-order 
Species Code 
Gradient A Gradient B 
Channa asiatica CHA* 6.3/0.1 8.3/0.1 
Siniperea roulei COR 12.5/0.3 33.3/0.3 
Ctenogobius sp. CTS 93.8/12.6 68.8/4.4 25.0/2.0 
Hypseleotris swinhonis HYS 12.5/0.2 9.3/0.1 
Mastacembelus aculeatus MAA 21.9/0.6 6.3/0.2 
Macropodus chinensis MAC* 16.7/«0.1 
Odontobutis obscurus ODO 12.5/0.9 37.5/2.7 91.7/7.9 
Siniperca chuatsi SIC 3.1/«0.1 25.0/1.2 75.0/2.4 
Siniperca obscura SIO 6.3/0.2 50.0/0.6 
Siluriformes 
Glyptothorax fukiensis GLF* 6.3/0.9 
Liobagrus styani LIS 31.3/1.0 31.3/0.6 
Mystus macropterus MYM* 8.3/0.1 
Pelteobagrus fulvidraco PEF 25.0/3.9 100.0/21.1 
Pseudobagrus albomarginatus PSA 3.1/«0.1 37.5/0.9 66.7/4.1 
Silurus asotus SIA 3.1/0.1 25.0/0.5 41.7/0.7 
Synbranchiformes 
Monopterus albus MOA 31.3/0.5 18.8/0.4 33.3/0.9 
Beloniformes 
Hyporhamphus intermedius HEI* 8.3/0.9 


*: Rare species occurring in less than two sites, not included in statistical analysis. 


Table2 Three-way ANOVA results of the spatial and temporal variations in fish species richness within Qingyi Stream 





Factors df SS MS F P Student-Newman-Keuls 

Stream orders 4 2.3 0.6 38.9 EE Ist<2nd=3rd<4th<5th 
Seasons 3 0.1 0.02 1.9 ns 

Gradients 1 0.1 0.1 4.0 v 

Ordersxseasons 12 0.2 0.01 0.9 ns 

Seasonsx gradients 3 0.02 0.01 0.5 ns A<B 
Orders*gradients 4 0.6 0.1 9.5 Er 

Ordersxseasonsxgradients 12 0.1 0.01 0.8 ns 


Factors are stream orders (first-fifth), seasons (January, April, July, October) and gradient (A, B; with and without the effect of a large dam, respectively). 


*: P<0.05, **: P<0.01, ns P>0.05. Student-Newman-Keuls test: P<0.05. 


Species richness significantly increased with stream 
orders, with the exception of no significant difference 
between second- and third-order streams. Gradient B 
had more species than that of Gradient A. A 
significant interaction effect was observed for stream 
order by gradient (Table 2). Within Gradient A, the 
lowest and highest species richness occurred in first- 
order and fifth-order streams, respectively, and species 
richness did not vary significantly among second-, 
third- and fourth-order streams. Within Gradient B, 
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however, the highest species richness was observed in 
fourth-order streams, while first-, second- and third- 
order streams were not significantly different in 
species richness (Figure 2). When comparing species 
richness within the same-order streams between the 
two gradients, significant differences were only 
observed in fourth-order streams (t-test, F=0.23, 
P«0.01), but not in first- (F=0.4, P=0.54), second- 
(F=1.5, P-0.21) or third-order streams (F=7.2, 
P=0.19) streams. 
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Figure 2 Fish species richness along two longitudinal gradients 
by stream orders within Qingyi Stream (upper: gradient A; 
lower: gradient B) 

Different lowercases represent significant differences in species richness 

(One-way ANOVA, P<0.05) . 


Fish assemblages 

The two-way crossed ANOSIM results suggested 
that fish assemblages significantly differed across stream 
orders but not seasons, and were uniform for both 
longitudinal gradients of A and B (Table 3). Within 
gradient B, significant differences in inter-order fish 
assemblages were almost exclusively observed, with the 
exception of fourth- vs. fifth-order. Whereas, fish assem- 
blages in fourth-order streams were similar with those in 


Table 3 R-values and their significance levels (P) for 
comparisons of fish assemblages among stream orders and 








seasons using ANOSIM 
Gradient A Gradient B 
R P R P 

Order 0.50 0.001 0.66 0.001 
Ist vs. 2nd" 0.06 0.440 0.60 0.029 
Ist vs. 3rd 0.38 0.012 0.42 0.036 
Ist vs. 4th 0.29 0.049 1.00 0.029 
Ist vs. 5th 1.00 0.001 0.99 0.001 
2nd vs. 3rd 0.06 0.420 0.80 0.029 
2nd vs. 4th” 0.06 0.370 1.00 0.029 
2nd vs. 5th 0.98 0.001 0.96 0.002 
3rd vs. 4th" 0.19 0.296 1.00 0.029 
3rd vs. 5th 0.90 0.001 0.93 0.002 
4th vs. Sth” 0.85 0.002 0.06 0.310 
Season -0.25 0.988 -0.10 0.971 


One-way ANOSIM was conducted for gradient (B) without the effect of 
Chencun Hydropower Station, two-way crossed ANOSIM was conducted 
for gradient (A) disrupted by the station and both gradients (Both). 
Significance level was accepted at P<0.05 (in bold). *: opposite significance 


level occurring between the two gradients. 
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second- and third-order but not fifth-order streams within 
gradient A (Table 3). For fish assemblages within the 
same-order streams, significant between-gradient differe- 
nces were observed in fourth-order streams, but not in 
first-, second, and third-order streams (Table 4). 


Table 4 R-values and their significance levels for compari- 
sons of fish assemblages among seasons and between two 
gradients disrupted by a large dam or not using ANOSIM 








R P 
Season Gradient Season Gradient 
Ist-order —0.33 —0.25 0.876 0.815 
2nd-order —0.04 0.00 0.562 0.617 
3rd-order —0.31 —0.73 0.800 0.988 
4th-order —0.38 0.75 0.790 0.020 
5th-order —0.19 0.872 


Two-way crossed ANOSIM was conducted for first- to fourth-order streams, 
one-way ANOSIM was conducted for fifth-order streams. Significance level 
was accepted at P«0.05 (in bold). 


Based on NMS analysis, fish assemblages within 
gradient B varied substantially from headwaters to main- 
stream, despite the overlap between fourth- and fifth- 
order streams (Figure 3). In contrast, within gradient A, 
substantial overlaps were observed between neighboring 
orders from first- to fourth-order streams, while fish 
assemblages in fourth-order streams were distinct from 
those in fifth-order streams. 
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Figure 3 Non-metric multi-dimensional scaling (NMS) by stre- 
am orders and dams for fish assemblage data in Qingyi Stream 
Black symbols, gradient disrupted by Chencun Hydropower Station; Open 
symbols, gradient without damming; Grey symbols, mutual part of both 
gradients. Normal triangle, reverse triangle, square, diamond, and rotundity, 


fish assemblages from first- to fifth-order streams, respectively. 


SIMPER analysis revealed that Z. platypus was 
more abundant in first-, second- and third-order streams 
and less in fourth- and fifth-order streams within gradient 
B. The occurrence and abundance of C. rarus, P. fasciata, 
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R. ocellatus, and S. nigripinnis determined how fish 
assemblages varied substantially from first- to fourth- 
order streams. Ctenogobius sp. and C. argus contri- 


buted most to the divergence in fish assemblages betw- 
een fourth- and fifth-order streams within Gradient A 
(Table 5). 


Table 5 Diagnostic species by SIMPER analysis for first- to fifth-order streams along gradients A and B within Qingyi Streams 











Gradient A Gradient B 
AvAbu. AvSim. Contri (%) AvAbu. AvSim. Contri (96) 

1st-order AvSim., 46.3% AvSim., 44.1% 

Zacco platypus 60.9 18.5 40.0 167.3 9.4 30.2 

Acrossocheilus fasciatus 39.6 10.2 22.1 17.3 1.5 17.9 

Vanmanenia stenosoma 16.3 5.1 11.0 

Cobitis rarus 23.5 24 19.3 
2nd-order AvSim., 51.6% AvSim., 52.9% 

Zacco platypus 63.1 16.8 32:5. 37.0 VAI 33.5 

Acrossocheilus fasciatus 15.3 9.7 18.9 15.0 3.1 24.7 

Ctenogobius sp. 14.4 7.4 14.4 

Parabotia fasciata 16.0 7.6 14.4 
3rd-order AvSim., 48.6% AvSim., 57.896 

Zacco platypus 44.8 11.9 24.6 95.8 15.9 27.6 

Rhodeus ocellatus 22.1 9.9 20.5 69.7 11.4 19.7 

Ctenogobius sp. 28.5 6.6 13.6 28.8 8.4 14.6 
4th-order AvSim., 51.0% AvSim., 64.3% 

Ctenogobius sp. 24.6 11.9 23.4 

Zacco platypus 46.1 8.5 16.8 

Acrossocheilus fasciatus 12.6 6.1 11.9 51.8 22 8.1 

Sarcocheilichthys nigripinnis 36.5 5.7 8.8 

Abbottina rivularis 20.3 4.8 7.5 
5th-order AvSim., 52.9% AvSim., 52.9% 

Pseudobagrus albomarginatus 58.0 7.0 15.1 58.0 7.0 15.1 

Channa argus 16.8 4.5 92 16.8 4.5 9.2 

Carassius auratus 15.8 3.4 7.3 15.8 34 7.3 


The first three species contributing most to the average similarity within each order are shown. AvSim., average similarity (%); AvAbu., average abundance; 


Contri., contribution. 


Relationships between species and environment 

The first and second axes of the CCA ordination 
accounted for 49.7% and 12.4% of total variance in 
species richness and abundance among sampling sites, 
and for 54.7% and 13.6% of variance in species data and 
environmental relationship, respectively. From CCA, 
gradient A separated fish assemblages based on elevation, 
substrate, and flow velocity on the right, and wet width, 
water depth, discharge, and conductivity on the left. The 
species related to this gradient were E. ilishaeformis, P. 
pekinensis, S. nigripinnis, H. labeo, C. rarus, C. argus, S. 
chuatsi, P. parva, C. auratus, A. chankaensis, R. 
ocellatus, C. sinensis, Ctenogobius sp., A. fasciatus, and 
Z. platypus (left to right). Gradient B was caused by 
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water temperature and dissolved oxygen on the top and 
pH on the bottom, along which P oxycephalus, P. 
fasciata, A. chankaensis, P. vaillanti, M. albus, G 
taeniellus, and S. argentatus were distributed from top to 
bottom (Figure 4). 


DISCUSSION 


Our results showed that along a longitudinal 
gradient without the direct effect of a large dam, fish 
species richness was highest in fourth-order streams and 
showed no significant difference among first- to third- 
order streams. Following the increase in stream size and 
habitat diversity from headwaters to downstream, fish 
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Figure 4 Canonical correspondence analysis (CCA) diagrams 
for environmental variables and fish assemblages in Qingyi 
Stream 
Black arrows indicate the most important factors and grey arrows show the 


other factors. Species codes as in Table 1. 


species richness generally increases downstream (Gorm- 
an & Karr, 1978; Matthews, 1986). Based on fish speci- 
mens collected from 89 shallow (first- to third-order) 
tributaries within the Qingyi watershed, Yan et al. (2011) 
discovered that local species richness was determined by 
local habitat variables (e.g., wetted width and water 
temperature) but not tributary spatial variables (e.g., 
stream order and link magnitude). Yan et al (2010) also 
revealed that stream order was not the optimal frame- 
work explaining spatial variation in fish assemblages 
within Puxi Stream, a tributary above Chencun Reservoir 
in the Qingyi watershed. The results in our study that fish 
species richness showed no substantial variation across 
first- to third-order streams appears to support fish 
species responding to local habitat features but not to 
stream classification schemes of humans (Matthews, 
1998). In addition, our findings that the highest species 
richness occurred in fourth-order, not fifth-order, streams 
is consistent with the river continuum concept (RCC) 
claiming that maximum species diversity often occurs in 
mid-sized, not large, streams. Oberdorff et al (1993) also 
found a decline in fish species richness in the lower 
reaches of the Seine River, which was explained by 
anthropogenic disturbances decreasing habitat diversity. 
In addition to species richness, fish also exhibit an 
alteration in zonation along the upstream-downstream 
gradient (Lasne et al, 2007; Miranda & Raborn, 2000). 
One of the most recognized theoretical classifications 
dividing running waters based on fish is longitudinal fish 
zonation proposed by Huet (1959); the brown trout, 
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grayling, barbel and bream zones exhibit alteration, in 
turn, from headwaters to lower reaches. However, asso- 
ciated with continuous, not abrupt, variation in the 
natural environments of many streams, fish assemblages 
vary by “addition” rather than by “replacement” (Matthe- 
ws, 1998), which suggests a deficiency in the application 
of this zonation concept within these streams. In our 
results, despite some species only being collected in one 
ordered stream (e.g., P. oxycephalus and C. molitorella in 
first-order streams, C. carpio, M. piceus, and H. 
leucisculus in fifth-order streams), most fish species 
were distributed in at least two different ordered streams. 
Based on our SIMPER analysis, along the longitudinal 
gradient without the effects of a large dam, changes in 
the relative abundances of Z. platypus, C. rarus, P. 
fasciata, R. ocellatus, and S. nigripinnis contributed most 
to longitudinal species replacement. Our results suggest 
that “addition”, not “replacement”, was the main underl- 
ying mechanism explaining species distributions from 
upstream to downstream in Qingyi Stream. This spatial 
variation in species distributions could be explained by 
the relationship between species and the local habitat 
environment based on CCA. The spatial distributions of 
fish species in this study were affected by a series of 
local habitat environmental factors, such as elevation, 
substrate, water depth, current discharge, wetted width, 
and conductivity. At the watershed scale, Yan et al (2011) 
reported that fish assemblage structure in this study area 
was related to the combined effects of local habitat (i.e., 
elevation, substrate, and water depth) and spatial stream 
position (i.e. magnitude link and confluence link), 
suggesting that fish species distribution was associated 
with their ecological requirements, such as habitat prefe- 
rence and trophic ecology. 

In Qingyi Stream, fish assemblages presented diffe- 
rent patterns in their spatial variations between the two 
longitudinal gradients, which could be viewed as the 
effect of Chencun Hydropower Station. Fish showed 
significantly lower species richness and different species 
composition within the gradient disrupted by the dam 
compared with the gradient without the large dam. 
However, the substantial variations in fish assemblages 
between gradients were only observed in fourth-order 
streams. Numerous studies have revealed that dams may 
impact stream fish by blocking migratory pathways 
(March et al, 2003), decreasing current velocity upstream 
(Bennett et al, 2002), altering natural flow regime 
(Bonner & Wilde, 2000), and decreasing water temp- 
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erature downstream (Clarkson & Childs, 2000). In this 
study, when considering first- to fourth-order streams 
within the two gradients independently, a total of 18 fish 
species (e.g., D. tumirostris, E. bambusa, S. hollandi) 
were collected in gradient A only, and eight species (e.g., 
O. bidens, H. bleekeri, S. parvus) were more frequent 
and abundant in gradient A than in gradient B. This may 
be explained by two causes resulting from the effects of 
large dams on fish assemblages. Firstly, most of the 18 
species only occurring in gradient A, including =H. 
bleekeri, E. ilishaeformis, H. maculatus, and S. roulei, 
inhabited the lower reaches of the streams. The Chencun 
Hydropower Station may play a role in constraining their 
upstream movements, resulting in the local extinction of 
these fish in areas above the station. Secondly, dams may 
habitat | conditions 


modify local upstream and 
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